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1^-^ ' The equilibrium properties of particles adsorption is investigated theoretically. The model relies on 

f*"^ I a free energy formulation which allows to generalize the Maxwell-Boltzmann description to solutions 

, for which the bulk volume fraction of potentially adsorbed particles is very high. As an application 

' we consider the equilibrium physical adsorption of neutral and charged particles from solution onto 

^ ' two parallel adsorbing surfaces. 

■ PACS numbers: 



^ • I. INTRODUCTION 

O 

,^ ^ The adsorption phenomena, due to the electrochemical interaction between the particles of a system and a surface, 

' is present in many experimental set up, such as the adsorption of a perfect gas on a surface [J or of charged particles 
in an electrolyte In many physical or chemical system, a better understanding of the theoretical equilibrium 



o 



and references therein) would thus be useful to interpret the 



o 



experiments. Several works have been devoted to this question and various models of particles distribution have been 
proposed. Among these, many assumed a Maxwell-Boltzman particle distribution (see 00 0). for instance, 
Q study the ionic adsorption on a surface due to some electrochemical forces in order to determine the surface density 
I— 1 1 of adsorbed charges versus the thickness of the sample. This work helps to understand the thickness dependence of 

the anisotropic part of the anchoring energy experimentaly observed 6] in a nematic liquid cristal 4^ 7] 8^ . 
^ I However, a limit can be made about the Maxwell-Boltzman distribution. Actually, this distribution can only 
. correctly describe the distribution properties in the dilute regime. But even in this regime, the density is usually large 
' at the surface itself, exept when the affinity of the particle for the surface is weak. To overcome the restriction to the 
T— I , first limit, the dilute case, we propose to apply a free energy formalism to the study of the equilibrium properties of 
■ neutral and charged particles adsorption onto two parallel adsorbing surfaces. Another advantage of the free energy 
formalism lays in the fact that it leads to the generalized Poisson-Boltzmann equation introduced in [Tol | which takes 
into account the finite size of the ions. In that paper the behavior of electrolytes solutions closed to a charged surface 
was studied. In our work, the surface is rather charged by the adsorption of one of the two charges present in the 
. system. Within our framework we obtain the electric potential distribution from the generalized Poisson-Boltzmann 
equation, the correct equations for the bulk particle distribution and the density of the particles on the surface with 
respect to the thickness d of the sample. For small thickness, the d dependence of the electric potential and of the 
'"^ ' chemical potential are determined and it is found that the surface density is proportional to thickness. Whereas in 
' the limit of large d the electric potential and the surface density are independent of the thickness. It is nevertheless 
O clear that results obtained with the phenomenological, coarse-grained free energy formalism to systems approaching 

I molecular dimension can only be trusted as far as general trends are concerned. 
^ I It is important to note that in the context when several kind of particles are present in the system, as it would 
be in an electrolyte, we find a particle distribution different from the Fermi-Dirac like distribution introduced in 0. 
Actually, the Fermi-Dirac distribution takes naturally into account the occupation of the adsorption sites. Yet it 
misses the mixing entropy contribution which is present in our formalism. As a consequence, we show that in the 
limiting case of a weak electrolyte, the results of the Poisson-Boltzmann approach are recovered by our formalism but 
not by the Fermi-Dirac distribution. 

The paper is organized as follows. In section 2 we introduce the mean field free energy formalism for neutral particle 
in a isotropic fluid limited by two adsorbing surfaces. In the same section the case of the adsorption competition 
between two neutral particles is studied. In section 3 we generalize the free energy formalism to the study of the 
ionic adsorption in a isotropic fluid limited by two adsorbing surfaces, already studied by means of the Fermi-Dirac 
distribution in 3] and 9]. 
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II. NEUTRAL PARTICLE ADSORPTION 



A. Theory on a lattice 



Consider N neutral particles in a slab of thickness d delimited by two surfaces of area S. We divide the slab into 
discrete cells of size (the size of the particle) and each cell is limited to a single particle occupation. We call Nt 
the number of sites in the bulk and Ng the number of surface adsorption sites. In thermodynamic equilibrium, rif, 
and Us are the number of particles in the bulk and at the surface respectively. The volume fraction is then — ni,/Ni, 
and the surface density (j)s = Us/Nf,. The conservation of the total number of particle N — Us + ni, is also written: 

<i> = 20,^+</.(l-^) (1) 

which is valid if the bulk volume fraction is uniform. If this fraction is not uniform (see section 3), relation (|Hl 
becomes: 

Q I Ad-2a)/2 

$ = 20, - + - / (f>{x)dx (2) 

a «J-(d-2a)/2 

where $ is the total volume fraction. 



B. Free energy formalism 

We will use the free energy formalism for the particles adsorption. It has been first introduced by D. Andelman et 
al. (for a review see pj|) to describe the kinetic adsorption of surfactant. This theoretical approach was successfully 
applied to the kinetic of non-ionic and ionic surfactant adsorption as well as to the kinetic of surfactant mixture 
adsorption. In this formalism, the two equations describing both the diffusive transport of surfactant molecules 
from the bulk solution to the interface and the kinetic of adsorption at the interface itself are derived from a single 
functional. The scope of the present paper is to apply the free energy formalism to the study of the equilibrium 
properties of particles adsorption. 

Following prj l we write the total free energy as a functional of the volume fraction in the bulk (j){x) and the density 
at the interface (jjg 

F(6) Ad-2a)/2 

-^ = 2fM+ I f[cl^{x)]dx (3) 

where the bulk free energy density is written: 



-(d-2a)/2 



/(</>) = _<;A:T[((/)ln0+(l-0)ln(l-0)]-^02_~^l 




and the surface free energy density is equal to : 

In </., + (!- cj)s) In (1 - (/.,)] - a^s ~ \<i>l ~ /I<^s| (5) 

The parameters a accounts for the energetic preference of the particle to absorb on the surface. /? is the lateral 
interaction between two adjacent charges. Note that the main difference with the free energy introduced in [nj in the 
context of surfactant adsorption lies in the presence of the exact entropic term in (0)) rather than on an approximate 
term. Actually, below as well as above the critical micellar concentration, the free chain surfactant solution is always 
dilute, so that a good approximation for the entropic term is 5 = 01n0 — 0. Another difference is due the effect on 
the finite volume: in Q and /I is the chemical potential at equilibrium. Its value is not imposed by an external 
reservoir but is determined by the conservation equation . In the infinite volume case considered in TTl , the chemical 
potential is imposed by an external reservoir localized at infinity. This last condition imposes the equilibrium bulk 
volume fraction. 
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The variation of F with respect to (/i(a;),that is 

SF 



54) {x) 

yields the bulk equilibrium volume fraction 



(6) 



(t){x)=(t)= ^- — ^ (7) 

in which we introduce the dimensional quantities /3 = fi/kT and /i = Jl/kT. At the surface the condition 

d(ps 



yields the equilibrium adsorption isotherm 

1 



1 + e-t^-'^+0<t>a 



with a = a/kT. 

Considering the case /? = 0, we find the Fermi-Dirac distribution that can also be written 



^ ?!)+(l-(/))e-" 

For /3 = 0, the number of particles conservation allows to compute analytically the chemical potential 



2$ 



(9) 



(10) 



e-^ = -a+Va2 + 6 (11) 

where 

[$ - (1 - 2a/d)] + ($ - 2a/d) e" 



(12) 



(1 - $) 

b=^ ^e" (13) 

$ 

Note that for d = 2a we can check that e^^ = e" (1 — $) /$ which leads to the expected result 0s — <i>, since all 
particles are localized on the surfaces. 



1. Case <E> >> 2a/d 

The condition 

1 > $ >> 2a/d (14) 

yields for the chemical potential 

- - (-) 

and a bulk volume fraction of : 

(/) = 4>. (16) 

Condition \1A\ implies a negligible variation of the equilibrium volume fraction after the adsorption process. At the 
surfaces, 

$ 

(17) 



^ $ + (l-$)e-" 

is independent of the size d. Even if the surface density of particles is large, the sample is large enough to ensure that 
the bulk volume fraction does not change. This result can also be obtained from the particle number conservation 
equation ([Ql $ = 20s| + cj){l — ^) which in the limit d — > oo gives </) = $. 
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2. Case $ << 2a/d 



This condition corresponds to a dilute regime <&<<!, where << 1 so that /x « Ini/). The surface coverage can 
then be written 

(18) 



or from the particle number conservation 

d$ — 2a4>s 



(19) 



d$~2a(/)^ + (d-2a)e-("+'3^=)/^ 
The dilute regime for /? = is of some interest. For a dilute solution the chemical potential (|lll) is approximatively 

_ 1 + 2a/d(e" - 1) 



leading to : 



^ l + 2a/fi(e"-l) 

and, at the surface 



2a/(i + e-" (1 - 2a/d) 



(20) 



(21) 



(22) 



Note that this relation can also be obtained from equation H19|l . 
If d << 2ae", is negligible and (|22|l becomes 

4.^^l (23) 

that is the surface coverage increases linearly with the size of the sample. Note that in paper ^ the same expression 
was written 

N 

(bs ~ d (24) 

where N is bulk density of particles in the absence of adsorption and Ns is the surface density of sites. With the 
identifications $ = iVa^ and Ns — the two expressions (|23|l and (|24|l coincide. But, as the authors of did not 
introduce a lattice, expression l|24|) leads to the unphysical result 0s ~* in the limit d ^ 0, since the correct limit 
d — > 2a is hidden. Note that this problem will be even more apparent in the case of the ionic adsorption. 
In the opposite limit d » 2ae" from (|21|l we deduce for the volume fraction in the bulk 

which imposes : ^s^^ << $. Actually, from we see 

^s^^ (26) 

from which we deduce 

2a , , 

(l)s^«^ (27) 

In this case the sample is large enough so that the volume fraction can be considered as constant even when the 
surface density is large. The system is then equivalent to an infinite system coupled to an external reservoir, this last 
one keeping the volume fraction constant. We see a crossover between a regime where the surface coverage increases 
linearly and another regime in which the surface coverage is independent of thickness. 
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C. Equilibrium distribution of two kinds of neutral particles 

In this section, we consider an infinite system composed of two neutral species which can both adsorb on a flat 

surface. 

We generalize the free energy formulation designed in the preceding section by writing the bulk contribution of the 
density free energy in kT units 

kT 

f{4>A,<t>B) = ^ {(7)Ahl0A + 0Bln<^B + (1 - 0^ - (^B)ln(l - - 0b) (28) 



~(t>A-^4>'B-^(l)A(l>B- ^J.A4>A- IJ.B(I>B^ (29) 



/3a ,2 /5b ,2 



kT 

fs{(t^s,A,4's,B) = ^ {(/)s,Aln0s,yl + ^i's.B ln(/)s,B + (1 - (/-s./l - </'s,B)ln(l- (30) 



where e is an interaction between the two species. At the surface, we have : 

(l>s,A - (t>s,B) - a^0s,A - aB(f>s,B - ^4>l,A " ^4'1,B (31) 
-£(t^s,A<f>s,B - f-lA(t>s,A - ^^b4's,b} (32) 

Note the presence of the mixing entropic term (1 — 0^ — (/)b) In (1 — — (/)b) in these two expressions. This terms 
is very important since it avoids that two particles of different kind sit at the same place in the lattice. Its absence 
would lead to the FD distribution. 

Minimizing the free energy, we obtain in the bulk; 

1 



and 



whereas at the surface we have: 



and 



1 -)- g-iHA+0A4>A+S<l>B) 

l-<t>A 



(33) 
(34) 



"'^ 1 -I- Q-{lJ.A+aA+l3A4>,,A+e4>,,B) ^ ' 



\ -|- g^-{fJ.B+ClB+0B4>a,B+e4>s,A) 

We thus see that the distributions of the two species are not independent of each other, due to the mixing entropy. 
Suppose now that a a » olb- We find 



'^'^'^ ~ et^A+CA 

and 

1 



« 1 (37) 



(38) 



showing that only one specie adsorbs, the other staying in the bulk. One can check that (j)s,A + 0s, s is always smaller 
than one. 

Now, let us compare our result with the Fermi-Dirac distribution. In such a context, the distribution for the two 
species are: 

1 



1>A,B 



1 + e~^-*'^ 
for the bulk and 

1 



,A,B 



(39) 



(40) 



\ -\- g— MA.B— o;a,b 

for the surface. The two distributions are now completely independent. In particular for the >> o-b the sum 
<i>s,A + 4's,B is not guarantee to be smaller than 1. This example shows the importance of taking the mixing entropy 
into account when more than one specie are present. 
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III. IONIC ADSORPTION 



The power of the free energy formaUsm can also be apphed to the ion distribution in an isotropic fluid hmitcd 
by two adsorbing surfaces. As explained in Q and 0, this system has already been considered by several authors. 
Actually, the ionic adsorption has been invoked to explain the thickness dependence of the anisotropic part of the 
anchoring energy of the interface between a substrate and a nematic liquid crystal. 

Consider a slab of thickness d with two identical adsorbing flat surfaces that adsorb only positive ions. Obviously 
the liquid is globally neutral. However, due to the selective ionic adsorption there is a distribution of charges yielding 
a non uniform electric potential V (x) across the sample. Since the surfaces are identical i.e. the affinity of the positive 
ions for the surfaces are identical, the potential is symmetric V (x) = V {—x) and E = —dV/dx is vanishing at the 
middle of the sample. 

The total free energy for a symmetric electrolyte in the mean held approximation has already been introduced in 
[lo| in the context of the adsorption of large ions from a solution of inflnite size to a charged surface. In our case it 
is rather the adsorption phenomena which charges the surfaces. Within the mean-fleld approximation, the total free 
energy in the bulk / = u~Ts can be written in terms of the local electrostatic potential in kT unit ^'{x) = eV (x) /kT 
and the ions volume fraction 0^ (x). The electrostatic energy contribution is 



kT 



u = 



dx 



r2 



dx 



(41) 



where e is the dielectric constant of the solution, are the equilibrium chemical potential of the two ions and 
Lb = is the intrinsic length of the problem. Note that we use the same system of unit as which is different 

EkTo?_ rpj^g term in the left hand side of H41I) is the self energy of the electric 



from the one of [lOl| where £r = 



held, the next two terms are the electrostatic energy of the ions. For the sake of simplicity we do not introduce 
additional steric interaction. The entropic contribution is 



Ts 



kT 



dx I 



.ln( 



Int 



+ (1-0+ In (1-0+ -(/._)] 



(42) 



The flrst two terms represent the tanslational entropy of the ions and the last term the entropy of mixing i.e. the 
entropy of the solvant molecules. 

At the interface itself, the total free energy is obtained by adding an electrostatic contribution (jf^tps to equation 



kT 

fs{<j>t.i^s) = ^{[0^1n0 

The total free energy of the system is then: 

^=2A(0+,V^, 



1- 



ln(l- 



(rf-2a)/2 



{x),ij{x)]dx 



-(d-2a)/2 

The ions number conservation imposes the equality between the two chemical potentials /i+ = 
of the bulk free energy with respect to (f>^ yields the volume fraction of the ions in the bulk: 



and 



(43) 
(44) 

fi^ ^ fi. The variation 
(45) 

(46) 



where 



h{ip,ii) = 1 + 2e^coshV' 



(47) 
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Note that the distribution in the bulk is very different from the FD distribution which reads 

The variation of the bulk free energy with respect to ip yields the modified Poisson Boltzman equation introduced in 
M : 



„9 , sinh , ,„x 



This equation is similar to the expression given in 9]. The difference lies in the formula for h{^jj,ii) which is in (9j : 
h{ilj,fi) = 1 + 2e^coshV' + e"^''. 

The variation of the total free energy with respect to ip {x = d/2) = ips yields the requirement of the overall charge 
neutrality: 



dx 



= T^-^. (50) 

x=d/2 



where Lg = o^ekT/e^ is a caracteristic length of the surface introduced in 0. Minimizing the surface free energy 0s 
yields 

(51) 



1 + e'^^"-'' 

which is a Fermi-Dirac distribution. 

Note that a FD distribution is obtained for when the electrostatic potential is very high » 1 since in this 
case: 

(52) 



whereas 

<t>^ - --"'<t>- ^ (53) 

It is instructive to analyze the behaviour of the system for semi infinite and very thin sample and then to compare 
the prediction of our model to the results obtained with the FD distribution. 



A. Infinite volume limite 



In the infinite volume limit d 



oo, 



we have 



(54) 



whereas 



^FD 



(55) 



due to the lack of the mixing entropy. The generalized PB equation can be solved numerically to find the electric 
field and the ions distributions across the sample. This has been done in ^3 where interesting curves can be found. 
It is only in the dilute case that the two chemical potential coincide and are equal to the Boltzman one : 



'Boltz 



(56) 



In this case, we obtain the surface potential 




(57) 
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and the surface coverage : 

1 



, -V^ (58) 
1 + 



in agreement with 



1. Small volume limit 
In the finite volume case we have two conservation laws 

d'^^+ciy„(,_2„)/2l + 2eMcMV') ^^^^ 

I Ad-2a)/2 n,+4, 

$ = - / (60) 

i-(<i-2a)/2 1 + 2e^cMV') 

Let us consider the small volume limit d 3a. Physically, we cannot consider a smaller bound since the negative 
charges are not adsorbed In particular, it is not possible to take the limit d ^ 0. 
The relation ()6Ufl becomes (assuming e'^'" ^ 1, which will be justified later): 

^ d-2a e^'+'^o 

1—-, TT- 61) 



which leads to 



Considering the approximation : 



e-'^'" (62) 



■0s ~ "00 - -r-0s (63) 
and the fact that is very small, we can assume that : 

i^s ~ ipo (64) 

in l(5T)) . so that : 



Plugging this result in H59|) and using the fact that e" ^ 1 yields the electric potential : 



^V-o pa/2 . / d_ 



1 _ ^ _ $ 



confirming our assumption e'^° ^ 1. 

For the chemical potential we readily obtain: 



1— 2a_(j)A; 2a_(j) 

d V 

and for the surface coverage we find : 



(65) 



(66) 



' ^ ^e-"/2. (67) 



(68) 
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These results are very different from the one obtained by Barbero et al. with the Fermi Dirac distribution which are 

~ — (69) 



$ (70) 

These two last relations do not lead to the result obtained with the Maxwell-Boltzman distribution used in the dilute 
regime. In other words, the Fermi- Dirac distribution in the limit of a small concentration does not lead to the correct 
Maxwell- Boltzman result. 

On the contratry, our equations H66|) and (|67|l . for $ <C 1, yields : 



and 



2a (d - 2a) 



$e~"/2 (72) 



that are the results obtained with a Maxwell-Boltzmann distribution. 

An estimation of the parameter was given in Jij for a typical nematic liquid crystal (e « 6 for an organic liquid) 
limited by two glasses. The adsorption energy was evaluated a « 6 and for a typical molecule of radius R ~ 40 A one 
has Lb ~ 30A. The surface density was found to be d dependent for thickness smaller than 300 A. 



IV. CONCLUSION 



In this paper, we have proposed a free energy formalism to describe the phenomenon of surface adsorption of 
neutral and charged particles as well. This free energy formalism has led to the equilibrium particles distribution for 
the case of physical adsorption of neutral and charged particles from solution onto two parallel adsorbing surfaces. 
In particular, we have found the correct equations for the electric potential and the equilibrium charge distribution 
with respect to the thickness of the electrolyte sample in case of high bulk concentration and we recover the results 
obtained with the Maxwell-Boltzmann distribution in the limit of small concentration. 

We are aware that our model relies on some strong assumptions, in particular the adsorbed particles are confined 
to a monomolecular layer whereas multilayer adsorption is frequently observed. Moreover, we have assumed that the 
surface is homogenous which is obviously not the case in general. 

Nevertheless, one of the advantages of the free energy formalism is that it relies on a minimization principle, avoiding 
in this way the introduction of ad-hoc distributions, and allowing straightforwardly the description of multi particles 
adsorption. Moreover it can be extended to take into account the particles interactions in the bulk and at the surface 
itself, that are usually not considered. This generalization needs further investigation. 
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